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封 结构 泄漏 流动 气 热 耦 合 研究 
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摘 要 : 本 文采 用 气动 传 热 耦合 计算 方法 研究 了 轮 载 封 严 冷 气 和 叶 冠 泄漏 流动 对 涡轮 性 能 以 及 热 负 荷 的 影响 。 研 究 结果 表明 
端 区 泄漏 流动 会 引起 流动 损失 并 导致 涡轮 做 功能 力 下 降 , 不 同 股 泄漏 流动 对 涡轮 效率 的 影响 可 以 近似 满足 线性 须 加 规律 。 固 
体 部 件 热传导 对 叶片 气动 负荷 几乎 没有 影响 ， 而 由 轮 载 封 严 冷气 造成 的 温度 梯度 会 显著 改变 壁面 换 热 特性 ， 此 时 必须 采用 气 
热 粮 合 计算 才能 准确 预测 热 锋 荷 。 
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Aerothermal Analysis of a Turbine with Endwall Leakage Flow 


AN Yu-Ge! LIU Huo-Xing! & ZOU Zheng-Ping 


(1.National Key Laboratory of Science and Technology on Aero-Engines Aero-Thermaldynamics, School of Energy and Power 
Engineering, Beihang University, Beijing 100191, China) 
Abstract: This study presents the results of the aerothermal analysis of an axial turbine with endwall leakage flow using conjugate 
heat transfer computation (CHT). As it turns out, a linear relationship of performance penalty between different flow systems is 
observed. The results of thermal analysis demonstrate that although heat conduction has little impact on the aerodynamic 


characteristics, it strongly affects the temperature distribution especially in the presence of coolant and large temperature gradients. 
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等 人 四 采用 非 定 常 的 实验 和 数值 方法 研究 了 封 严 气 
体 对 于 涡轮 性 能 以 及 二 次 流动 带 来 的 影响 ， 并 且 分 

航空 发 动机 遵循 布雷 顿 循 环 ， 提 高 涡轮 前 温度 析 研 究 了 冷气 预 旋 角度 对 冷却 效果 以 及 附加 的 流动 
从 而 增 大 循环 功 并 提高 热效率 是 研制 高 推 重 比 高 性 损失 的 作用 ，Ameri 和 Bunker9g 通 过 绝热 壁 和 等 温 
能 先进 航空 发 动机 的 重要 途径 。 现 代 航 空 发 动机 的 壁 的 CFD 计算 获得 转子 叶 顶 的 对 流 换 热 系数 分 布 ， 
涡轮 燃气 温度 已 经 远 远 超过 了 部 件 材 料 的 耐 温 能 指出 了 叶 尖 热 负 蓓 与 泄漏 流动 结构 密切 相关 ， 精 确 
力 ， 因 此 在 航空 发 动机 中 广泛 采用 冷却 技术 保证 涡 把 握 流动 结构 特征 是 准确 预测 热 负 蓓 的 前 提 ; Haller 
轮 部 件 的 安全 工作 ， 尤 其 是 在 转 静 子 之 间 的 盘 腔 会 和 Hilditcht "研究 了 叶 冠 外 表面 的 对 流 换 热 特性 。 以 
采用 封 严 冷气 来 防止 燃气 入 侵 趾 ， 并 在 叶 顶 采用 带 上 热 分 析 方 法 仅仅 通过 CFD 计算 获得 绝热 壁 温 或 
冠 结 构 减 少 泄漏 流量 。 据 统计 在 发 动机 中 引 自 压气 流体 域 边界 换 热 特性 ， 如 果 需 要 得 到 实际 部 件 温度 
机 的 冷却 气体 高 达 主 流 流量 的 20% 中 ， 其 中 仅 涡 轮 还 需要 进行 额外 的 热传导 计算 ， 因 此 被 称 作 非 耦合 
封 严 冷气 就 可 以 消耗 整 机 耗 油 率 的 6%P1, 由 于 叶 冠 热 分 析 方 法 。 随 着 计算 技术 的 进步 ， 气 动 传 热 耦合 
泄漏 流体 不 做 功 也 会 引起 效率 随 泄漏 流量 的 线性 下 计算 方法 (Conjugate heat transfer，CHT) 也 逐步 在 
降 外 。 先 进 航 空 发 动机 对 涡轮 性 能 越 来 越 高 的 要 求 涡轮 叶片 热 分 析 外 ， 优 化 内 冷 系统 设计 外 上 获得 应 
迫使 研究 者 们 开始 关注 泄漏 流动 对 涡轮 性 能 和 热 负 用 。 这 些 研究 加 深 了 我 们 对 涡轮 端 区 密封 结构 泄漏 
人 答 的 影响 ， 采 用 精细 化 设计 减 小 冷气 量 和 附加 的 气 流动 对 主流 流 场 和 热 负 蓓 的 影响 的 认识 。 然 而 在 实 
动 损失 。 际 工作 中 ， 涡 轮 主流 流动 会 同时 受到 轮 载 叶 冠 等 多 

随 着 技术 的 进步 越 来 越 多 的 研究 者 针对 泄漏 流 股 泄漏 流动 的 影响 ， 在 封 严 冷气 引起 的 强 温 度 梯度 


动 对 涡轮 流 场 和 热 负荷 的 影响 展开 研究 ， 例 如 Ong 的 影响 下 固体 域 热传导 现象 也 会 影响 到 部 件 的 换 热 
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特性 ， 随 着 涡轮 负荷 的 进一步 增 大 ， 主 流 和 泄漏 流 
动 系统 之 间 的 相互 作用 、 流 场 和 温度 场 之 间 的 耦合 
影响 必 将 更 加 强烈 。 这 就 要 求 我 们 对 传统 研究 针对 
单一 的 轮 载 或 叶 冠 泄漏 流动 进行 的 分 析 得 到 的 结论 
以 及 采用 单纯 CFD 计算 获得 涡轮 部 件 热 负 荷 分 布 
规律 的 适用 性 和 可 推广 性 进行 审视 。 欧 盟 的 第 六 框 
架 计 划 中 的 MAGPID 和 AITEBU9 项 目 也 开始 展开 
类 似 的 探讨 ， 针 对 计 及 真实 涡轮 环境 的 端 区 泄漏 流 
动 引 起 的 气动 传 热 相关 现象 进行 研究 。 
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气 入 侵 。 同 时 对 于 本 涡轮 转子 盘 的 后 表面 直接 暴露 
于 燃气 中 ， 因 此 在 封 严 冷气 和 主流 燃气 的 影响 下 在 
转子 盘 内 部 形成 了 非常 强 的 温度 梯度 ， 流 场 和 温度 
场 之 间 也 会 有 强 相 互 作用 。 在 计算 中 主流 进口 给 定 
总 温 总 压 ， 封 严 冷气 进口 给 定 流量 和 温度 ， 在 流体 
域 转动 部 件 和 静止 部 件 之 间 采 用 挫 混 面 模型 ， 流 体 
和 固体 域 边 界 保证 温度 和 热流 连续 ， 在 固体 域 的 外 
边界 给 定 温度 分 布 。 计 算 网 格 如 图 2 所 示 ， 主 流 道 采 
用 分 块 结构 化 网 格 ， 由 16 个 网 格 块 构成 ， 轮 坑 、 叶 


本 文 研究 工作 采用 气动 传 热 耦 合 的 计算 方法 有 
完了 轮 载 封 严 冷气 和 叶 冠 泄漏 流动 对 某 单 级 涡轮 流 
场 和 性 能 的 影响 以 及 不 同 流动 系统 之 间 流 动 特征 的 
非 线 性 破 加 效应 ;针对 轮 载 封 严 冷气 对 部 件 热 负 荷 
影响 进行 了 研究 ， 并 在 此 基础 上 分 析 了 强 温度 梯 
度 条 件 下 的 涡轮 部 件 温 度 场 与 流 场 的 耦合 作用 机 制 
及 其 对 传统 非 耘 合 热 分 析 方 法 计算 精度 带 来 的 影 
响 。 

1 数值 模拟 方法 
1.1 计算 程序 

本 文 工 作 采用 商用 程序 CFX 进 行 气动 传 热 耦合 
计算 ， 流 体 域 求解 雷诺 平均 的 NS 方程 ， 采 用 SST 
(shear stress transporb) 消 流 模型 对 方程 进行 封 朵 , 使 
用 7-Re, 转 损 模 型 预测 转换 并 考虑 其 对 流动 和 换 热 
带 来 的 影响 ， 在 固体 域 通过 求解 能 量 方程 来 模拟 热 
传导 。 对 于 非 耦合 热 分 析 方 法 , 首先 采用 绝热 壁 CFD 
计算 得 到 参考 温度 7 的 分 布 ， 并 通过 式 (1) 算 得 壁面 
温度 四 分布 ， 其 中 67 为 可 以 使 壁面 产生 热流 的 微小 
温 升 ， 再 采用 给 定 壁 面 温度 的 方式 进行 第 二 次 CFD 
计算 获得 热流 密度 g 的 分 布 , 进而 可 以 通过 式 (2) 得 到 
流体 域 边界 的 对 流 换 热 系数 ， 随 后 以 采用 第 三 类 边 
界 条 件 以 绝热 壁 温 为 参考 温度 ， 并 给 定 边界 换 热 系 
数 来 求解 固体 域 导热 方程 ， 即 可 获得 涡轮 部 件 温度 
分 布 。 


一 


(1) 


(2) 
1.2 计算 模型 及 网 格 

采用 带 有 盘 腔 及 叶 冠 结构 的 菜 发 动机 单 级 跨 音 
涡轮 为 算 例 进 行 研究 ， 该 涡轮 结构 如 图 1 所 示 。 涡 轮 
设计 转速 17000r/min, 设计 点 膨胀 比 为 1.8, 载荷 系数 
为 1.95。 冷却 气体 通过 迷宫 封 严 进入 转 静 子 之 间 的 盘 
腔 ， 降 低 涡轮 盘 热 负荷 ， 并 最 终 注 入 主流 道 防 止 燃 


Cl 


冠 封 严 腔 和 固体 域 采 用 了 非 结构 化 的 六 面体 网 格 ， 
在 近 壁 面 进行 了 加 密 ， 保 证 y+ 小 于 1， 三 维 计算 网 格 
总 单元 数 为 830 万 。 在 计算 分 析 之 前 进行 了 网 格 无 关 
性 验证 ， 结 果 表 明 当 继续 增加 网 格 数目 时 涡轮 输出 
功 和 叶 冠 泄漏 流量 及 叶片 平均 温度 几乎 无 变化 。 


图 1 涡轮 结构 示意 图 
Fig. 1 Turbine model geometry 


图 2 计算 网 格 


Fig. 2 Computational grid 
2 计算 结果 及 讨论 


2.1 泄漏 流动 对 性 能 的 影响 
为 分 析 轮 载 和 叶 冠 泄漏 流动 对 涡轮 性 能 的 影响 
及 各 个 泄漏 流动 系统 之 间 的 相互 作用 ， 本 文 分 别 对 
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理想 主流 道 流 动 、 叶 冠 泄漏 影响 下 的 流动 、 轮 载 封 
严 冷 气 影响 下 的 流动 以 及 叶 冠 泄漏 和 轮 载 封 严 同时 
作用 下 的 流动 进行 了 计算 研究 。 图 3 给 出 了 各 种 不 同 
几何 复杂 度 计算 获得 的 涡轮 效率 ， 定 量 地 比较 了 轮 
载 叶 冠 泄漏 流 动 对 涡轮 性 能 造成 的 有 影响。 值得 注意 
的 是 ， 在 气 冷 涡轮 发 展 过 程 中 众多 的 学 者 和 研究 机 
构 以 及 企业 基于 多 种 掺 混 做 功 模 型 对 涡轮 效率 的 定 
义 进 行 了 研究 ， 然 而 就 目前 而 言 尚 未 对 效率 的 定义 
方式 达成 一 致 站， 本 文 的 分 析 研 究 采 用 了 工业 上 广 
泛 采 用 的 效率 定义 方式 ， 如 式 G) 所 示 ， 分 母 为 燃气 
和 冷气 的 等 灶 膨 胀 功 ， 分 子 为 扭矩 乘 以 转速 即 涡轮 
实际 输出 功 。 从 图 中 可 以 看 出 轮 载 封 严 冷 气 会 导致 
涡轮 效率 下 降 1.18%， 而 叶 冠 泄漏 流动 会 引起 高 达 
2.29% 的 效率 下 降 。 当 同时 存在 着 叶 冠 泄漏 和 轮 载 封 
严 时 效率 的 下 降 与 两 者 分 别 造成 的 效率 损失 之 和 非 
常 接近 ， 仅 有 0.02% 的 效率 差别 。 这 说 明 在 本 算 例 轮 
载 和 叶 冠 泄漏 流动 系统 对 性 能 的 影响 可 以 看 做 是 线 
性 又 加 的 。 同 时 需要 进一步 指出 的 是 对 于 本 文 计算 
模型 ， 叶 冠 泄漏 流量 高 达 主 流 流量 的 2%， 按 照 传统 
的 叶 冠 损失 理论 ， 泄 漏 流量 造成 的 做 功能 力 损失 及 
效率 下 降 等 于 其 流量 百分比 ， 因 此 叶 冠 泄漏 的 做 功 
损失 是 本 算 例 泄漏 流动 最 大 的 损失 来 源 。 


EN TIw 
PF mCpeATss + mCpAT,, G) 
叶 冠 泄漏 + 轮 载 封 严 做 功能 力 损 失 
Amn=-0.02% An=-2.04% 

叶 冠 泄漏 流动 An=-2.29% 

轮 载 封 严 An=-1.18% 

主流 道 计算 n=91.98% 

一 
86 87 88 89 90 91 92 


效率 
图 3 泄漏 流动 对 涡轮 效率 的 影响 

Fig. 3 Influence of leakage flow on efficiency 
2. 2 泄漏 流动 对 热 负 荷 的 影响 

端 区 泄漏 流动 在 对 主流 燃气 流动 结构 以 及 涡轮 

性 能 造成 显著 影响 的 同时 ， 也 影响 着 涡轮 部 件 热 负 
荷 的 分 布 。 如 图 4 即 为 通过 气 热 耦合 计算 获得 的 涡 
轮 固体 部 件 的 无 量 纲 温度 分 布 ， 无 量 纲 温度 9 可 以 
通过 式 (4) 求 出 。 由 于 封 严 冷气 和 主流 燃气 之 间 存 在 
着 温差 ， 导 致 涡轮 盘 内 部 形成 了 很 强 的 温度 梯度 ， 
尤其 是 存在 于 转子 叶 根 的 温度 梯度 会 对 结构 强度 产 
生 不 利 影响 。 在 这 种 情况 下 涡轮 盘 的 热传导 现象 也 
会 显著 地 影响 封 严 气体 的 流动 换 热 特性 ， 而 换 热 特 
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生 的 改变 又 会 反 过 来 影响 温度 分 布 ， 从 而 形成 一 个 
遇 合 问题 。 图 5 给 出 了 采用 气 热 耦 合计 算 和 绝热 壁 
计算 获得 的 转子 叶片 5%、50% 和 95% 叶 高 三 个 截面 
的 等 炉 马赫 数 分 布 ， 等 恼 马赫 数 可 以 通过 式 (5) 计 算 
获得 ， 其 中 pz 是 转子 进口 的 平均 总 压 。 计 算 结 果 
表明 绝热 壁 计算 获得 的 气动 负荷 分 布 与 通过 气 热 耦 
合计 算 获得 的 结果 几乎 完全 相同 ， 这 说 明 叶 片 热 传 
导 对 气动 负荷 的 影响 非常 小 ， 在 气动 设计 与 计算 中 
完全 可 以 将 之 忽略 。 
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图 4 涡轮 部 件 温度 分 布 
Fig. 4 Temperature distribution of the turbine 
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一 上 一 Adiabatic Wall 95% Spa 
一 上 CHT 50% span 

一 一 Adiabatic Wall 50% Spa 
一 一 CHT 5% span 

一 -Adiabatic Wall 5% Spanw 
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到 5 转子 叶片 表面 等 炉 马 赫 数 分 布 


Fig. $ Isentropic Mach number distribution 
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传统 的 非 耦合 热 分 析 方 法 通过 给 定 壁面 温度 获 
得 表面 对 流 换 热 系数 或 努 赛 尔 数 等 参数 的 分 布 来 评 
价 涡轮 部 件 热 负 荷 ， 必 须 通过 额外 的 热传导 计算 才 
可 以 获得 温度 分 布 ， 这 种 计算 方法 无 法 考虑 真实 回 
体 壁 面 温度 对 于 流体 换 热 特性 的 影响 。 图 6 对 比 了 
采用 CHT 计算 、 非 耦合 计算 以 及 绝热 壁 计算 获得 的 
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转子 盘 盘 腔 一 侧 温 度 径 向 分 布 。 可 以 看 出 在 实际 工 
作 中 涡轮 盘 的 温度 远 高 于 绝热 壁 温 。 由 于 涡轮 盘 的 
后 表面 直接 暴露 在 高 温 燃 气 中 ， 热 传导 现象 导致 了 
盘 腔 内 壁 温度 升 高 ， 进 而 改变 了 壁面 换 热 参数 。 图 
7 给 出 了 CHT 计算 和 非 耦 合计 算 获 得 的 无 量 纲 化 的 
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燃气 温度 ,Ty 和 工分 别 为 壁面 两 侧 第 一 层 单元 流体 
和 固体 温度 分 布 ，L4ny 和 Ln, 分别 为 网 格 单元 宽度 。 
在 轴 向 温差 驱动 下 平板 发 生 热传导 ， 热 流 密度 g 可 
以 通过 公式 (6) 计 算 获得 , 联 立 公式 (2) 和 公式 (6) 即 可 
得 到 壁面 温度 Tv 的 表达 式 如 式 (7)。 正 是 由 于 非 耘 


对 流 换 热 系数 分 布 ， 可 以 看 出 由 于 非 耦 合计 算 相 当 
于 获得 的 是 参考 温度 附近 的 换 热 特性 ， 因 此 较 实 际 
值 发 生 了 偏差 。 
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图 6 转子 盘 温 度 分 布 


Fig. 6 Temperature distribution on the disk 
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Normalized heat transfer coefficient 
图 7 转子 盘 对 流 换 热 系数 分 布 

Fig. 7 Heat transfer coefficient of the disk 
图 8 详细 给 出 了 采用 CHT 计 算 方 法 和 非 耦 合计 
算 获 得 转动 部 件 的 温度 分 布 之 比 ， 在 这 里 可 以 认为 
气 热 厅 合计 算 的 结果 更 接近 物理 实际 。 结 果 表 明 非 
耦合 计算 方法 高 估 了 盘 腔 表面 的 温度 分 布 ， 并 直接 
导致 了 转子 叶 根 的 温度 预测 值 偏 高 ， 而 叶 根 处 的 热 
负 蓓 恰恰 对 机 械 强 度 分 析 非 常 重要 。 为 进一步 分 析 
存在 较 大 轴 向 温度 梯度 条 件 下 的 转子 盘 气 动 传 热 厅 
合 问题 ， 将 其 简化 为 一 个 厚度 为 L 的 平板 ， 如 图 
9 所 示 ， 其 中 ,为 盘 腔 壁面 温度 ，7, 为 转子 盘 后 壁 


合计 算 获 得 了 偏 低 的 壁面 对 流 换 热 ， 因 此 导致 了 对 
壁面 温度 预测 偏 高 。 
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图 8 气 热 耦合 计算 与 非 耦 合计 算 结果 对 比 


Fig. 8 Temperature ration distribution 
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图 9 热传导 对 转子 盘 温 度 的 影响 
Fig. 9 Conjugate heat transfer of the disk 
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3 结 论 


本 文采 用 了 气动 传 热 耦 合计 算 方法 研究 了 轮 载 
和 叶 冠 泄漏 流动 对 涡轮 性 能 和 热 负荷 的 影响 ， 并 比 
较 分 析 了 涡轮 部 件 的 耦合 和 非 耦 合 热 分 析 方 法 对 热 
负荷 预测 的 准确 度 。 主 要 结论 如 下 : 

1) 轮 载 和 叶 冠 泄漏 流动 会 引起 流动 损失 并 导致 
涡轮 做 功能 力 下 降 ， 不 同 复杂 度 的 泄漏 流动 系统 对 
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hh 


生 能 的 影响 符合 线性 县 加 规律 。 

2) 固体 部 件 热传导 对 热 负 荷 分 布 具有 重要 影 

响 ， 尤 其 是 在 轮 载 封 严 冷气 造成 的 温度 梯度 的 作用 
下 这 一 现象 更 为 突出 ， 然 而 热传导 对 叶片 气动 负荷 
几乎 没有 影响 ， 因 此 对 于 气动 性 能 工程 设计 和 计算 
采用 绝热 壁面 假设 是 合理 的 。 
3) 当 涡 轮 部 件 内 部 温度 梯度 较 大 时 ,热传导 会 
影响 壁面 温度 进而 影响 对 流 换 热 系 数 ， 必 须 采 用 气 
动 传 热 耦 合计 算 才 能 获得 符合 物理 实际 的 温度 分 
布 。 
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